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1 Introduction 

In metallic clusters the electron wave function is confined in all three spacial 
dimensions, resulting in a discrete set of energy levels. In quantum wires the 
confinement of the electron wave function is limited to two spacial dimensions, 
which gives rise to the formation of a set of one-dimensional energy bands, 
which we refer to as modes, or channels. 

The influence of these quantum modes on the properties of metallic quan- 
tum wires has been studied experimentally by forming point contacts between 
metal electrodes, using scanning tunnelling microscopes, mechanically con- 
trollable break junctions, or related techniques. Where the primary obser- 
vation of quantum size effects in clusters was based on cluster abundance 
spectra, for metal point contacts there is a somewhat analogous statistical 
method, which consists in recording a histogram of conductance values ob- 
served for large numbers of contacts. For simple, monovalent metals these 
histograms show a number of pronounced peaks close to multiples of the 
quantum of conductance, 2e 2 /h. 

The interpretation of the histograms is less straightforward compared to 
that of the cluster abundance spectra, since the former are the result of a 
combination of atomic structural features, and of the quantisation of the 
electron states in ID subbands. Various experimental methods have recently 
been used to investigate the role of quantised modes in atomic-size contacts. 
These include measurements of the mechanical force on the contacts, the use 
of characteristic features in the current-voltage relation of superconducting 
contacts, measurements of the thermopower, of the voltage-dependence of the 
conductance and measurements of the shot noise intensity. With the help of 
these methods it can be shown that the conductance of atomic size contact 
for simple metals (Au, Na, etc.) is indeed carried by a well defined set of 
quantum modes. However, in general the conductance is determined by the 
number and character of the valence orbitals of the metal atoms forming the 
contact. 

A very direct analogy between the physics of metal clusters and the 
physics of metallic point contacts becomes visible when recording conduc- 
tance histograms for the alkali metals and for contacts larger than just a few 
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atoms. A large number of peaks is observed, which have the same origin as 
the "magic numbers" in cluster abundance spectra. The observations suggest 
that the electronic quantum mode structure influences the mechanical sta- 
bility of the nanowire, giving preference to those diameters which correspond 
to a filled "shell" of conductance modes. 

The outline of this chapter is as follows. First a brief introduction is given 
of the natural formalism for discussing electron transport in ballistic conduc- 
tors: the Landauer theory. After introducing the experimental techniques, 
which are used for studying ballistic point contacts in metals, the experimen- 
tal observations for the conductance of atomic-scale contacts are presented. In 
order to obtain a full description in terms of the quantum modes for conduc- 
tance, several recently developed techniques are reviewed, which go beyond 
straightforward measurement of the conductance. A brief discussion is given 
of an unusual atomic geometry for gold contacts, which evolve into a chain 
of freely suspended atoms. Then shell filling effects in sodium nanowires are 
discussed in the context of the influence of the conductance modes on the 
total energy of the system. The chapter ends with an outlook on promising 
new developments. 



2 Landauer Theory for Ballistic Conductance 



The metallic point contacts and nanowires, which we will consider, are all 
of atomic size, much smaller than all characteristic scattering lengths of the 
system. In particular, the electron mean free path for elastic scattering on 
defects and impurities near the contact is assumed to be much longer than 
the contact size. The only elastic scattering considered is the scattering by 
the walls forming the boundary of the system. Also, it is assumed that the 
probability for scattering events, which change the spin and phase of the 
electron wave function, is negligible. Following the standard approach (for 
reviews, see |l],g,[|) we will assume that the system can be schematically 
represented as in Fig. ffl. The connection between the ballistic system and 
the outside world is represented by electron reservoirs on each side of the 
contact, which are held at a potential difference eV by an external voltage 
source. When the leads connecting the reservoirs to the contact are straight 
wires of constant width, there is a well defined number of conducting modes 
in each of these wires, say N and M for the left and right lead, respectively. 
In a free electron gas model the modes are simply plane waves, which can 
propagate in the current direction (to the left and right) and are standing 
waves in the perpendicular directions. The modes can be labelled by an index 
corresponding to the number of nodes in the perpendicular direction. The 
numbers N and M arc limited by the requirement that the energy of the 
modes is lower than the Fermi energy. 
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Fig. 1. Schematic representation of a ballistic point contact. The reservoirs on the 
left and right fully absorb the incoming electron waves. The lead on the left has 
a width, which admits N conductance channels, and the one on the right has M 
channels 



The conductance of the system can now be simply expressed as |l],|]|| , 

G = ^Tr(ttt) , (1) 

where e is the electron charge, h is Planck's constant, and t is an NxM matrix 
with matrix element t mn giving the probability amplitude for an electron 
wave in mode n on the left to be transmitted into mode m on the right 
of the contact. It can be shown that the product matrix tH can always be 
diagonalised by going over to a new basis, consisting of linear combinations of 
the original modes in the leads. Further, the number, N c , of non-zero diagonal 
elements is only determined by the number of modes at the narrowest cross 
section of the conductor ||,|| . Equation (Q) thus simplifies to 

2e 2 ^ 
G = X^>' W 

n=l 

where T n — \t nn \ and the index refers to the new basis. 

Under favourable circumstances, all transmission probabilities, T„, can 
be close to unity. As an example, for a smooth ( "adiabatic" ) and long wire 
the modes at the narrowest cross section couple exclusively to a single mode 
at either side of the contact, and the expression for the conductance further 
simplifies to 

G = N C G , (3) 

where Go = 2e 2 /h is the conductance quantum. For increasing diameter of 
the contact the number N c increases each time a new mode fits into the 
narrowest cross section. This number is limited by the requirement that 
the kinetic energy for motion in the perpendicular direction is smaller than 
the Fermi energy. For a 2-dimensional (2D) system this can be expressed as 
(H 2 /2m)(TrN c /W) 2 < Ef, with W the width of the contact, which leads to 
iV c = Int(2PF/AF), with Af the Fermi wavelength. For a 3D metallic contact 
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Fig. 2. Schematic top and side view of the mounting of a MCB, where the metal to 
be studied has the form of a notched wire (1), which is fixed onto an insulated elastic 
substrate (3) with two drops of epoxy adhesive (4) very close to either side of the 
notch. The substrate is mounted in a three-point bending configuration between the 
top of a stacked piezo element (5) and two fixed counter supports (2). This setup is 
mounted inside a vacuum can and cooled down to liquid helium temperatures. Then 
the substrate is bent by moving the piezo element forward. The bending causes the 
top surface of the substrate to expand and the wire to break at the notch. Typical 
sizes are L ~ 20 mm and it ~ 0.1 mm 

N c rs (7ri?/Ap) 2 , with R the contact radius. This quantisation of the conduc- 
tance was first observed in experiments on 2D electron gas devices by van 
Wees et al. |j and by Wharam et al. 0, where Af — 400 A is much larger 
than the atomic scale. 



3 Experimental Techniques 

The experimental tools for fabricating atomic-scale contacts are mostly based 
on a piezoelectric actuator for the adjustment of the contact size between two 
metal electrodes. Standard Scanning Tunnelling Microscopes (STM) are of- 
ten used for this purpose P,p|,^0[ . The tip of the STM is driven into the 
surface and the conductance is recorded while gradually breaking the con- 
tact by retracting the tip. The first experiment of this type was reported by 
Gimzewski and Moller pT| . 

A practical tool for the purpose of studying metallic quantum point con- 
tacts is the Mechanically Controllable Break-junction (MCB) technique fll2[ |. 
The principle is illustrated in Fig. 0. By breaking the metal, two clean frac- 
ture surfaces are exposed, which remain clean due to the cryo-pumping action 
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of the low-temperature vacuum can. This method circumvents the problem 
of surface contamination of tip and sample in STM experiments, where a 
UHV chamber with surface preparation and analysis facilities are required 
to obtain similar conditions. The fracture surfaces can be brought back into 
contact by relaxing the force on the elastic substrate, while a piezoelectric 
element is used for fine control. The roughness of the fracture surfaces results 
in a first contact at one point, and experiments usually give no evidence of 
multiple contacts. In addition to a clean surface, a second advantage of the 
method is the stability of the two electrodes with respect to each other. From 
the noise in the current in the tunnelling regime one obtains an estimate of 
the vibration amplitude of the vacuum distance, which is typically less than 
10~ 3 A. The stability results from the reduction of the mechanical loop which 
connects one contact side to the other, from centimetres, in the case of an 
STM scanner, to ~ 0.1 mm in the MCB. 

Conductance properties of atomic-size contacts can even be studied by 
still simpler methods. When breaking the contact between two regular wires 
under ambient conditions [|13| or switching a commercial relay |14f and when 
recording the evolution of the conductance with a time resolution of the order 
of microseconds, the conductance can be observed to decrease by atomic steps 
during the last stages of the contact breaking. Such methods allow rapid 
accumulation of histograms of conductance values for a large number of such 
scans. 



4 Conductance of Atomic-Scale Contacts 

Figure g shows some examples of the conductance measured during break- 
ing of a gold contact at low temperatures, using an MCB device. The con- 
ductance decreases by sudden jumps, separated by plateaux, which have a 
negative slope, the higher conductance the steeper. Some of the plateaux are 
remarkably close to multiples of the conductance quantum, Go ; in particular 
the last plateau before loosing contact is nearly flat and very close to 1 Go- 
Closer inspection, however, shows that many plateaux cannot be identified 
with integer multiples of the quantum unit, and the structure of the steps is 
different for each new recording. Also, the height of the steps is of the order 
of the quantum unit, but can vary by more than a factor of 2, where both 
smaller and larger steps are found. Drawing a figure such as Fig. @, with grid 
lines at multiples of Go, guides the eye to the coincidences and may convey 
that the origin of the steps is in quantisation of the conductance. However, in 
evaluating the graphs, one should be aware that a plateau cannot be farther 
away than one half from an integer value, and that a more objective analysis 
is required. Still, it is clear that we can use these fairly simple techniques 
to produce and study atomic-scale conductors, for which the conductance is 
dominated by quantum effects. The interpretation of graphs as in Fig. will 
be the subject of this and the following sections. 
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The fact that the atomic structure and orbital character of the electron 
modes is important for an interpretation of the conductance can already be 
deduced by comparing the conductance curves for various metals. Figure 
shows three examples of conductance curves for aluminium atomic-size con- 
tacts, and Fig. H shows results for platinum. In both cases, in particular for 
aluminium, we find that many plateaux have an anomalous slope: the con- 
ductance increases when pulling the contact, in contrast to the results for 
gold. For aluminium, the last plateau before breaking is still close to one unit 
of conductance, but one frequently observes the conductance diving below 
this value, and then recovering to nearly 1 Go, before contact is lost. Plat- 
inum, on the other hand, has a last conductance value, which is usually of 
order two times larger and also the size of the jumps is somewhat larger. 

For all metals the transition between the plateaux is very sudden and 
sharp. Such sudden transitions would not be expected in a model which 
describes the conductance in terms of a set of conduction modes, which are 
gradually pinched-off by reducing the contact diameter. Indeed, the jumps 
find their origin in sudden rearrangements of the atomic structure of the 
contact fll7|Jl8| . Upon stretching of the contact, the stress accumulates elastic 
energy in the atomic bonds over the length of a plateau. This energy is 
suddenly released in a transition to a new atomic configuration, which will 
typically have a smaller contact size. Such atomic-scale mechanical processes 
were first described by Sutton and Pethica Q and by Landman et al. [M . 




100 150 200 250 300 



Piezo-voltage (V) 



Fig. 3. Three typical recordings of the conductance G measured in atomic size 
contacts for gold at helium temperatures, using the MCB technique. The electrodes 
are pulled apart by increasing the piezo-voltage. The corresponding displacement 
is about 0.1 nm per 25 V. After each recording the electrodes are pushed firmly 
together, and each trace has new structure. (After J.M. Krans ]15|) 
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Fig. 4. Three examples of the conductance measured in atomic size contacts for an 
aluminium MCB junction at 4.2 K, as a function of the piezo-voltage. (After Eaj) 
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Fig. 5. The conductance for a platinum junction at 1.3 K as a function of piezo- 
voltage for two successive scans. (After [p.6[) 



The first direct proof for atomic rearrangements at conductance steps 
was provided in an experiment by Rubio, Agrait and Vieira pjj ], where the 
conductance for atomic-size gold contacts was measured simultaneously with 
the force on the contacts (Fig. ||). The stress accumulation on the plateaux 
and the coincidence of the stress relief events with the jumps in the conduc- 
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tance can be clearly distinguished. Presently, the experiment has only been 
reported for gold at room temperature. At low temperatures, the evidence 
for the atomic structure related nature of the jumps comes from analysis of 
the dynamic behaviour of the jumps. Generally, hysteresis is observed in the 
position of the jumps when retracing the curve immediately after a jump is 
found (Fig. R). When increasing the bath temperature or the current through 
the contact (which indirectly heats the contact) the width of the hysteresis is 
seen to gradually decrease until it is reduced to zero, and spontaneous jumps 
between the two conductance values are observed, which have a thermally 
activated behaviour |18| , [22p3|| . These observations find a natural interpreta- 
tion in terms of jumps between two distinct atomic configurations for the 
contact, separated by an energy barrier. Various recent molecular dynamics 
simulations confirm this scenario [HH,G3,Efl.E7| . 



5 Histograms of Conductance Values 

The atomic configuration of a contact adjusts itself in response to the ex- 
ternally applied stress, and evolves depending on the starting configuration 
of the contact at larger size. The fact that each conductance curve differs in 
many details from previous curves reflects the fact that the atomic configu- 
ration of the contact for given conductance is different in each run. However, 
as observed in the previous section, there appears to be a certain preference 
for conductance values near integer multiples of the quantum unit. For gold 
near 1 Go this is immediately obvious from the examples in Fig. |3|. A general 
and objective method of analysis was introduced [g£|[29J, which consists in 
recording histograms of conductance values encountered in a large number 
of runs. Figure |3| shows a histogram for gold measured using a room temper- 
ature STM under UHV conditions |3JJ. Up to four peaks are found centred 
near the first four multiples of Go- For sodium in a low temperature exper- 
iment using the MCB technique, a histogram with peaks near 1, 3, 5 and 6 
times Go was observed [gj|. Figure || shows a similar result for potassium. 
The fact that peaks near 2 and 4 Go are absent points at an interpretation in 
terms of a smooth, near-perfect cylindrical symmetry of the sodium contacts. 
Sodium indeed forms a very good approximation to a free electron system, 
and the weakly bound s-electrons strongly reduce surface corrugation. This 
is also why the best abundance spectra for clusters, with pronounced features 
at magic numbers have been obtained for the alkali metals |3l| . For a model 
smooth, cylindrically symmetric contact with continuously adjustable con- 
tact diameter [ p2|]33] ], the conductance increases from zero to 1 Go as soon 
as the diameter is large enough, so that the first conductance mode is oc- 
cupied. When increasing the diameter further, the conductance increases by 
two units because the second and third modes are degenerate. The modes are 
described by Bessel functions (assuming a hard wall boundary potential) and 
the first mode is given by the in — Bessel function, which is not degenerate. 
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Fig. 6. Simultaneous measurement of force and conductance on atom scale point 
contacts for Au. The sample is mounted on a cantilever beam and the force between 
tip and sample is measured by the deflection of the beam using an Atomic Force 
Microscope (AFM). The measurements are done in air at room temperature. (From 



@D 



The second and third modes are the degenerate m = ±1 modes, followed by 
m = ±2 for further increasing contact diameter. The next mode which will 
be occupied corresponds to the second zero of the m — Bessel function, 
and is again not degenerate. Thus the conductance for such a contact should 
increase by 1, 2, 2 and 1 units, producing just the series of conductance values 
observed in the sodium experiment. The slight shift of the peaks in Fig. M 
below the integer values can be attributed to an effective series resistance 
due to back-scattering on defects near the contact. Model simulations [B4J of 
the histogram are in close agreement with the experiment, including the shift 
of the peaks. The fact that gold histograms show peaks at all the first four 



10 



J.M. van Ruitenbeek 



0.3 A 



3.5 



3.0 



2.5 



__'. 


T=2.2 K 


1 1 

\ 1 


■ 




82 83 84 85 86 87 88 89 90 91 



v p [V] 

Fig. 7. Measurements of the conductance for an atomic-size gold contact while 
sweeping the piezo-voltage forward and back over a single step. The curves have 
been recorded at three different temperatures while maintaining the sweep centred 
around the same jump. Clear hysteresis of the order of 0.1 A is observed at 2.2 K. 
At 3.1 K the hysteresis is half as large, and at 4.2 K it has disappeared. (After J.M. 
Krans pi) 



quantum values may be explained by a stronger deviation from cylindrical 
symmetry in gold, which then lifts the degeneracy of the modes. Model sim- 
ulations are indeed able to reproduce the shape of the gold histograms pfj . 



Not all metals show such pronounced histogram peaks near integer con- 
ductance values. The most clear-cut results are obtained only for monovalent 
metals. The alkali metals Li, Na and K show a histogram structure as repre- 
sented by Fig. 0, although for Li the shift of the peaks is somewhat stronger 
[ p8| . The noble metals Cu, Ag and Au show histograms as in Fig. f| US . The 
details such as the shift in position, the width and relative height of the peaks 
can be different depending on the experimental conditions p^ , p(j| , p9|j3^j3f^ , |37t . 
Most other metals only show a rather broad first peak, which reflects the con- 
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Fig. 8. Histogram representing the relative weight which each conductance value 
has in the experiments. The histogram is constructed from 227 conductance curves 
recorded while breaking Au contacts, using an STM under UHV at room temper- 
ature. (From @) 
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Fig. 9. Histogram of conductance values, constructed from G(V p )-curves measured 
for potassium at 4.2 K with an MCB device, involving several thousand individual 
measurements. The characteristic sequence of peaks (G = 1,3,5,6) is regarded as 
a signature for conductance quantisation. (From |Bq|) 



ductance of a single atom contact (see below). This peak can generally not 
be identified with an integer value of the conductance; for example niobium 
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shows a wide peak centred near 2.5 Go (Fig. noh and similar results have 
been obtained for Pb J39fl. On the other hand there are a few examples of 
multivalent metals, which show pronounced peaks in the histograms, among 
which aluminium Gol . As we shall discuss below, the histogram for Al throws 
doubt upon a straightforward interpretation of the histogram peaks in terms 
of conductance quantisation. For some systems other than simple metals ev- 
idence has been obtained for preferred conductance values near multiples of 
the quantum of conductance (carbon nanotubes, semimctals, metal oxides, 
etc.). These will not be discussed here, since the nature of the contact and the 
mechanism of contact formation is expected to be very different from that in 
ordinary metals. 

When we assume that the contact breaking process produces any effec- 
tive contact diameter with equal probability, then the histograms represent 
a derivative of conductance with respect to the effective diameter of the con- 
tact. It is instructive to calculate the integral of the histogram, as was first 
done by Gai et al. |3q]. Fig. O shows such a curve, obtained from a gold 
histogram similar to the one shown in Fig. 0. This curve is to be compared 
to conductance traces obtained for 2D semiconductor devices Mm , for which 
the width of the contact can directly and continuously be adjusted by the 
gate electrostatic potential. Compared to the latter, the conductance steps 
in Fig. O are poorly defined, with the exception of the first conductance 
quantum. Moreover, we will show that the first quantum feature results from 




G (2e7h) 

Fig. 10. Histogram constructed from 2400 individual conductance curves for a nio- 
bium sample. Each curve was recorded while stretching the contact to break, using 
the MCB technique at a temperature of 10 K, which is just above the supercon- 
ducting transition temperature. The conductance was measured using a DC voltage 
bias of 20 mV. (From B. Ludoph H) 
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the fact that our assumption mentioned above is not valid. The effective 
diameters produced during contact breaking are strongly influenced by the 
possible atomic configurations. As we will argue below, the step at 1 Go, 
corresponding to the strong peak in the histogram for gold, results from the 
formation of a chain of gold atoms during the last stages of contact breaking. 
Disregarding the first step in Fig. 11, we find that the conductance is not 



strictly quantised, as the probability of finding a contact with a conductance 
of, e.g., 2 Go is only twice that of finding 1.5 Go- However, the conductance 
is still determined by the quantum states, as described in Sect. |2] and is car- 
ried by a limited number of modes. We will show below how the quantum 
nature for monovalent metals is revealed by a tendency for the modes to open 
one-by-one as the contact becomes larger. 



6 The Character of the Conductance Modes Through 
a Single Atom 

Instead of discussing the average properties of many contacts, we now con- 
centrate on the simplest configuration, a single atom bridging the electrodes 
at either side, and consider the question of what the conductance for this 
atomic configuration will be. Ideally, we would like to know the number of 
modes contributing to the conductance and the transmission probability, T n 
for each of these. From a measurement of the conductance alone we can- 
not obtain this information, since the conductance gives only the sum of 
T n . Scheer ct al. Ell have introduced a method which allows us to obtain 
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this information from experiment. The method exploits the non-linearities in 
the current-voltage characteristic for contacts in the superconducting state. 
Other techniques, which give more limited information on the contribution 
of the various modes, will be discussed at the end of this section. 



6.1 Subgap structure in superconducting contacts 

The principle of the method introduced by Scheer et al. can be illustrated 
by considering first a contact having a single mode with low transmission 
probability, T<1. For T <C 1 we have essentially a tunnel junction, and the 
current -voltage characteristic for a superconducting tunnel junction is known 
to directly reflect the gap, A, in the density of states for the superconductor 
| )42| . As illustrated in Fig. |l2| a no current flows until the applied voltage 
exceeds 2A/e (where the factor 2 results from the fact that we have identical 
superconductors on both sides of the junction), after which the current jumps 
to approximately the normal-state resistance line. For eV > 2 A single quasi- 
particles can be transferred from the occupied states at Ep — A on the low 
voltage side of the junction to empty states at Ep + A at the other side. For 
eV < 2A this process is blocked, since there are no states available in the 
gap. 

However, when we consider higher order tunnel processes a small current 
can still be obtained. Figure |l2](b) illustrates a process, which is allowed for 
eV > A and consists of the simultaneous tunnelling of two quasiparticles 
from the low bias side to form a Cooperpair on the other side of the junc- 
tion. The onset of this process causes a step in the current at half the gap 
value, V — 2A/2e. The height of the current step is smaller than the step at 
2A/e by a factor T, since the probability for two particles to tunnel is T 2 . In 
general, one can construct similar processes of order n, involving the simulta- 
neous transfer of n particles, which give rise to a current onset at eV = 2A/n 
with a step height proportional to T". An example for n = 3 is illustrated in 
Fig. y_2p. This mechanism is known as multiple particle tunnelling and was 
first described by Schrieffer and Wilkins |43J. It is now understood that this 
is the weak coupling limit of a mechanism which is referred to as multiple An- 
dreev reflection Q[||[4|,[|^Q|9) . The theory could only be tested recently, 
since it requires the fabrication of a tunnel junction having a single tunnelling 
mode with a well-defined tunnelling probability T. For atomic size niobium 
tunnel junctions the theory was shown to give a very good agreement |5(j| . pl| . 
describing up to three current steps, including the curvature and the slopes, 
while the only adjustable parameter is the tunnel probability, which follows 
directly from the normal state resistance. 

Since the theory has now been developed to all orders in T Q||,|^|Q , 
Scheer et al. H] realised that this mechanism offers the possibility of ex- 
tracting the transmission probabilities for contacts with a finite number of 
channels contributing to the current, and is ideally suited to analysing atomic 
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Fig. 12. Illustration of the Multiple Particle Tunnelling (MPT), or Multiple An- 
dreev Reflection (MAR) processes. In each of the three diagrams the available 
quasi-particle states as a function of energy (vertical axis) at both sides of the 
tunnel barrier are given in the semiconductor representation. Dark grey are the 
occupied states, light grey the unoccupied states and the line in the middle of the 
gap represents the Cooper pair energy, which is separated by an energy A from 
the occupied and the unoccupied states. Applying an external electrical potential 
V across the junction shifts the states on the left side of the junction up by an en- 
ergy eV with respect to those in the right electrode. The ordinary superconducting 
tunnelling process is given in (a), which shows that a voltage V > 2A/e is required 
for single quasi-particles to cross the junction. The probability for tunnelling of a 
particle, determined by the transparency of the barrier, is T. This gives rise to the 
familiar jump in the current at eV — 2A in the current -voltage characteristic of a 
superconducting tunnel junction. When we consider higher order processes, the next 
order is represented in (b). This can be described as two quasi-particles crossing 
simultaneously to form a Cooper pair in the right electrode (MPT). Alternatively, 
the process can be regarded as an electron-like quasiparticle falling onto the barrier, 
which is reflected as a hole-like quasiparticle, forming a Cooper pair on the right 
(MAR). The two descriptions are equivalent, and give rise to a current step in the 
current-voltage characteristic of the junction at eV — 2A/2. The probability for 
the process is T 2 , since it requires the crossing of two particles, (c) shows the third 
order process, which involves breaking up a Cooper pair on the left, combining it 
with a quasiparticle, to form a Cooper pair and a quasiparticle on the right. It is 
allowed for eV > 2/1/3 and has a probability T 3 



size contacts. Roughly speaking, the current steps at eV = 2A/n are propor- 
tional to ^2 T^, with m the channel index, and when we can resolve sufficient 
details in the current-voltage characteristics, we can fit many independent 
sums of powers of T m 's. When the T m 's are not small compared to 1, all pro- 
cesses to all orders need to be included for a description of the experimental 
curves. In practice, the full expression for the current-voltage characteristic 
for a single channel from theory p6| , f47| , [48U49| is numerically evaluated for 
a given transmission probability T m (Fig. |13|, inset), and a number of such 
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curves are added independently, where the T m 's are used as fitting param- 
eters. Scheer et al. tested their approach first for aluminium contacts Ej . 
As shown in Fig. O, all Current-voltage curves for small contacts can be 
very well described by the theory. However, the most important finding was 
that at the last "plateau" in the conductance, just before the breaking of the 
contact, typically three channels with different T's are required for a good 
description, while the total conductance for such contacts is of order 1 Go, 
and would in principle require only a single conductance channel. Contacts 
at the verge of breaking are expected to consist of a single atom, and this 
atom would then admit three conductance channels, but each of the three 
would only be partially open, adding up to a conductance close to 1 Go- This 
very much contradicts a simple picture of quantised conductance in atomic 
size contacts, and poses the question as to what determines the number of 
channels through a single atom. 



6.2 Valence-Orbital-Based Description of the Conductance 
Modes 

Cuevas, Levy Yeyati and Martfn-Rodero H constructed a model to explain 
these results, using a tight binding calculation, for a geometry of two atomic 
pyramids touching at the apex through a single atom. They argue that it 
is very important to make the tight binding calculation self-consistent, by 
which they mean that local charge neutrality is maintained at each atomic 
site, by iteration and adjustment of the site energy for each individual atom 
in the model configuration. Figure [14] shows the results of their calculations 
of the density of states and the transmission probability for the various chan- 
nels as a function of energy. They find that the conductance channels can be 
described in terms of the atomic valence orbitals. Aluminium has a configu- 
ration [Ne]3s 2 3p 1 , and a total of four orbitals would be available for current 
transport: one s orbital and three p orbitals, p x ,Py and p z . They identify 
in their calculation three contributions, one which originates from a combi- 
nation of s and p z orbitals (where the z coordinate is taken in the current 
direction), and two smaller identical contributions labelled p x and p y . The 
degeneracy of these two channels is due to the symmetry of the problem, and 
can be lifted by changing the local environment for the central atom. The 
fourth possible channel, an antisymmetric combination of s and p z , is found 
to have a negligible transmission probability. Thus, their calculation confirms 
the experimental observation by Scheer et al. that three channels contribute 
to the conductance for a single aluminium atom. It was also found that the 
total conductance for the three channels is of order 1 Go . The results are very 
robust against changes in the atomic configuration; only the total conduc- 
tance varies somewhat between different choices for the atomic geometry. The 
qualitative features of the model agree with ab initio calculations for single 
atom contacts and a simplified structure for the environment p2,p3 54 59] . 
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Fig. 13. Current-voltage characteristics for four atom-size contacts of aluminium 
using a lithographically fabricated mechanically controllable break junction at 
30 mK (symbols). The right inset shows the typical variation of the conductance, or 
total transmission T — G/Go , as a function of the displacement of the electrodes, 
while pulling, and is similar to the curves for aluminium shown in Fig. 0. The 
bar indicates the approximate length scale. The data in the main panel have been 
recorded by stopping the elongation at the last stages of the contact (a-c) or just 
after the jump to the tunnelling regime (d) and then measuring the current while 
slowly sweeping the bias voltage. The current and voltage are plotted in reduced 
units, el jGA and eV/A, where G is the normal state conductance for each contact 
and A is the measured superconducting gap, A/e — (182.5 ± 2.0)nV. The left inset 
shows the current-voltage characteristics obtained from first-principles theory for 
a single channel junction [MjTj,kq,W9|] with different values for the transmission 
probability T (from bottom to top: T=0.1, 0.4, 0.7, 0.9, 0.99, 1). The full curves 
in the main panel have been obtained by adding several theoretical curves and 
optimising the set of T values. The curves are obtained with: (a) three channels, 
Ti=0.997, T 2 =0.46, T 3 =0.29 with a total transmission J2 T ™ =1-747, (b) two chan- 
nels, Ti=0.74, T2=0.11, with a total transmission "^2,T n =0.85, (c) three channels, 
Ti=0.46, T 2 =0.35, T 3 =0.07 with a total transmission J]T„ =0.88. (d) In the tun- 
nelling range a single channel is sufficient, here ^ T n = 21=0.025. (From Scheer et 
al. {§) 



The analysis, both theoretical and experimental, was extended to other 
metals Rpfl , by which it was shown that the number of conductance channels 
for an atom of a given metallic element depends on the number of valence 
orbitals. Figure [l5] shows conductance curves, similar to those in Fig. ||, for 
Pb, Al, Nb and Au, where at each point in the figure current-voltage curves 
as in Fig. 13 were recorded and fitted in order to determine the number of 



channels involved. This number of channels is indicated along the curves in 
Fig. [TBI The number is constant over a plateau in the conductance, where 
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Fig. 14. Localized orbital model for electrical conduction through one-atom con- 
tacts. The atomic valence levels (vertical bars above the figures) develop into bulk 
conduction bands (a) and (b) for Al and Pb, respectively. The panels (c) and (d) 
depict schematically the local density of states (LDOS) in eV -1 at the central atom 
of the model geometry for a single atom contact. The global energy dependence of 
the transmission coefficients T n is shown in (e) and (f). The dotted lines indicate 
the position of the Fermi level. The sp z mode is the best transmitted for both 
materials. The p x and p y modes are degenerate due to the symmetry of the model 
geometry. (From 



the transmission probability for every mode changes gradually. At the steps 
in the conductance the number of channels involved is usually found to jump 
to a smaller number. In tunnelling range, when the contact is broken and the 
distance is larger than 0.2 nm, the current-voltage characteristics can in all 
cases be described by a single channel, with a transmission probability which 
is given by the tunnelling resistance. The number of channels found for the 
smallest contacts, just before the jump to tunnelling is 1 for Au, 3 for Al and 
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Pb, and 5 for Nb. Note that gold is not a superconductor, and a special device 
was fabricated which allowed the use of proximity induced superconductivity 
[ p6| . The device is a nanofabricated version of a break junction, having a 
thick superconducting aluminium layer forming a bridge with a gap of about 
100 nm. This small gap was closed by a thin gold film in intimate contact 
with the aluminium. Superconducting properties were thereby induced in the 
gold film, and by breaking the gold film and adjusting an atomic size contact, 
the same subgap analysis could be performed. Both the Al and Au junctions 
were measured at temperatures of 100 mK, far below the superconducting 
transition temperatures. Pb and Nb were measured at 1.5 K. 

The number of channels and the total conductance at the last plateau 
before breaking found in the experiment agree very well with the theory, for 
which the results can be summarised in the following way. Single atom con- 
tacts for monovalent metals, including the noble metals Cu, Ag and Au, and 
the alkali metals, have a single valence orbital available for current trans- 
port, giving rise to a single channel with a transmission probability close to 
unity. The total conductance for such contacts is thus expected to be close 
to 1 Go- The experiment for gold shows indeed a single channel, but the to- 
tal conductance is a factor 2-3 smaller than predicted. This was tentatively 
attributed to the strong scattering in the nano-fabricated device, and other 
techniques, which will be discussed below, have found more close agreement 
with the predicted total conductance of 1 Go- For s-p metals, including Al 
and Pb, three channels give a noticeable contribution to the current. The 
fourth channel, i.e., the antisymmetric combination of s and p z is always 
nearly closed. The total transmission depends on the number of valence elec- 
trons. Pb has 4 valence electrons while aluminium has only three, and this 
has roughly the effect that the Fermi energy in Fig. U4 is shifted up to the 
peak in the p x ,Py channel for Pb. The total conductance for Pb is found to 
be nearly 3 Go in the calculations, while the conductance at the last plateau 
for lead does indeed start close to three quantum units, dropping to lower 
values upon further stretching of the contact. Niobium is a d-metal with a 
configuration [Kr]4d 5 s having 6 valence orbitals: 1 s and 5 d. The the- 
ory again predicts one combination with a negligible contribution and that 
the remaining five channels should add up to a total conductance of about 
2.8 Go, again in good agreement with the experiment. The overall agreement 
is good, but some variation is observed between contacts and the information 
is obtained on a limited number of contacts. Below techniques are discussed 
which give less detailed information, but directly show the average properties 
of the contacts, and confirm the picture presented here. 



6.3 Further Experimental Techniques 

Three other methods have recently been introduced in order to investigate the 
number of channels contributing to the conductance of atomic size contacts: 
the measurement of shot noise, conductance fluctuations and thermopower. 
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Fig. 15. Conductance curves measured as a function of contact elongation for Al, 
Pb, Nb and Au. The number of channels contributing to the conductance was 
determined at each point in the curves by recording the current-voltage relation 
and fitting the curves with the theory for superconducting subgap structure. The 
numbers along the curves in the figure indicate the number of channels obtained in 
this way. The number is constant over a plateau, and usually jumps to a smaller 
value at the steps in the conductance. (From |5q|) 



Shot noise is the result of the discrete character of the current due to the 
passage of individual electrons. It was originally found in vacuum diodes, and 
first discussed by Schottky in 1918 (5?|]. The passage of individual electrons 
can be regarded as a delta function of the current with time. The total current 
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is the sum of a random distribution of such delta functions, giving a time 
averaged current /, and a frequency spectrum of fluctuations which is white 
(up to very high frequencies), with a noise power equal to 2el. This shot 
noise can be observed, e.g., in tunnel junctions. 

For a perfect ballistic point contact, in the absence of back-scattering, 
i.e., all channel transmission probabilities are either 1 or 0, the shot noise 
is expected to vanish p^ , ^ , p0||6l| , |6^ , |63[ . This can be understood from the 
wave nature of the electrons, for which the wave function extends from the 
left bank to the right bank of the contact without interruption. When the 
state on the left is occupied for an incoming electron, it is occupied on the 
right as well and there are no fluctuations in this occupation number. In other 
words, the incoming electron is not given the choice of being transmitted or 
not, it is always transmitted when it enters an open mode. In order to have 
noise, the electron must be given the choice of being reflected at the contact. 
This will be the case when the transmission probability is smaller than 1 
and larger than 0. In single-channel quantum point contacts, shot noise is 
predicted to be suppressed by a factor proportional to T(\ — T), where T is 
the transmission probability of the conductance channel |pgp9 , pOj6 1| fi2| . This 
quantum suppression has recently been observed in point contact devices in 



a 2-dimensional electron gas [164 65 1. For a general multichannel contact the 



shot noise power is predicted to be 



P I = 2eVG o y^T n {l-T n ). (4) 



Since this depends on the sum over the second power of the transmission coef- 
ficients, this quantity is independent of the conductance, G — Go^2T n , and 
simultaneous measurement of these two quantities should give information 
about the channel distribution. 

When measuring shot noise on atomic size point contacts ]66|] it is nec- 
essary to work at low temperatures in order to reduce thermal noise, and to 
shield the contact carefully from external mechanical and acoustic vibrations. 
Using two sets of preamplifiers in parallel and measuring the cross-correlation 
of the noise for the two signals eliminates the noise of the preamplifier. Ex- 
perimental results for shot noise in gold point contacts for a number of con- 



ductance values are shown in Fig. 16. From the two measured parameters, G 
and Pi, one can determine at most two independent transmission probabili- 
ties. Instead, the results in Fig. Qq are compared to models which assume a 
certain evolution of the values for T n as a function of the total conductance, 
G. The full curve shows the simplest model, where all T n are either or 1, 
except for a single partially open channel, so that the total conductance is 
given by G/Go = N — 1 + Tn- The figure shows that the experimental data 
closely follow this behaviour, in particular at low conductance. In order to 
estimate the contribution of additional partially open channels the broken 
curves show the expected behaviour for an evolution of the channel trans- 
missions as illustrated in the inset. It shows that the deviation from a model 
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Fig. 16. Measured shot noise values for 27 gold contacts at 4.2 K with a bias current 
of 0.9 nA. Comparison is made with calculations in the case of one single partially 
transmitted mode (full curve) and for various amounts of contributions of other 
modes according to the model described in the inset (dashed curves). In the limit of 
zero conductance, these curves all converge to full shot noise, i.e., 2.9xl0 -25 A 2 /Hz. 
(From @) 



with the channels opening one-by-one is only 10% between 1 and 2 Go, and 
about 20% between 2 and 3 Go- The scatter in the data points is larger than 
the experimental error because each point was measured on a different con- 
tact, which can have a very different set of transmission values, and a purely 
systematic behaviour is not expected. Although we cannot determine each 
individual transmission value, we obtain information from the property that 
the noise increases the more channels are partially transmitted. The mini- 
mum possible noise for a given conductance value is represented by the full 
curve in Fig. |lfj, corresponding to a single partially open channel. The broken 
curves give an impression of the other partially open channels. 

There is a very strong suppression, down to 2% of the full shot noise 
value, for G — 1 Go- It needs to be stressed that this holds for gold contacts. 
There is a fundamental distinction between this monovalent metal and the 
multivalent metal aluminium, which shows no systematic suppression of the 
shot noise at multiples of the conductance quantum, and the intensities lie 
between 0.3 and 0.6 times 2el for G close to Go |6q]. 



Conductance Fluctuations. Interference between electron trajectories scat- 
tering on defects near atomic size metallic contacts, gives rise to dominant 
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contributions to the second derivative of the current with respect to bias 
voltage V, i.e., in dG/dV. This effect has the same origin as the well-known 
Universal Conductance Fluctuations (UCF) in diffusive mesoscopic conduc- 
tors pTfl . It was studied for point contacts an order of magnitude larger than 



the atomic size by Holweg et al. [p8| 



and by Ralph et al. |70[, and the the- 



ory was developed by Kozub et al. 71]. In experiments on gold contacts J72[ 
it was found that this voltage dependence is suppressed near multiples of the 
quantum value of conductance, n(2e 2 /h). By applying a constant modulation 
voltage at frequency 10 and measuring the current with lock-in amplifiers si- 
multaneously at uj and the second harmonic 2u the conductance and its 
derivative can be obtained during conductance scans as in Fig. |3j. From the 
combined data sets of many such curves one can construct a conductance 
histogram together with the average properties of dG/dV . It was found that 
dG/dV for a given conductance value has a bell-shaped distribution centred 
around zero and that the width of this distribution has a systematic variation 



with conductance. Figure 17 shows the standard deviation of the derivative 
of the conductance with bias voltage <jgv = ((dG/dU) 2 ), obtained from 3500 
curves for gold, where the conductance and the derivative of the conductance 
were measured simultaneously as a function of contact elongation. The con- 
ductance histogram for the same set of data is shown in the lower panel. 



The data for oqv display pronounced minima for G near multiples of Go . 
A similar effect can be recognised in the numerical calculations of Maslov et 
al. 173] , The explanation for this quantum suppression of the conductance 
fluctuations as presented in [|72| is illustrated in the inset of Fig. [L7|. The 
contact is modelled by a ballistic central part, which can be described by a 
set of transmission values for the conductance modes, sandwiched between 
diffusive banks, where electrons are scattered by defects characterised by an 
elastic scattering length l e . An electron wave of a given mode falling onto the 
contact is transmitted with probability amplitude t and part of this wave is 
reflected back to the contact by the diffusive medium, into the same mode, 
with probability amplitude a„ « 1. This back-scattered wave is then reflected 
again at the contact with probability amplitude r n , where T n = \t n \ = 
1 — |r n | 2 . The latter wave interferes with the original transmitted wave. This 
interference depends on the phase difference between the two waves, and this 
phase difference depends on the phase accumulated by the wave during the 
passage through the diffusive medium. The probability amplitude a n is a 
sum over all trajectories of scattering, and the phase for such a trajectory of 
total length L is simply kL, where k is the wave vector of the electron. The 
wave vector can be influenced by increasing the voltage over the contact, 
thus launching the electrons into the other electrode with a higher speed. 
The interference of the waves changes as we change the bias voltage, and 
therefore the total transmission probability, or the conductance, changes as a 
function of V. This describes the dominant contributions to the conductance 
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Fig. 17. (a) Standard deviation of the voltage dependence of the conductance ver- 
sus conductance for 3500 curves for gold. All data points in the set were sorted as 
a function of the conductance after which the rms value of dG/dV was calculated 
from a fixed number of successive points. The circles are the averages for 300 points, 
and the squares for 2500 points. The solid and dashed curves depict the calculated 
behaviour for a single partially-open channel and a random distribution over two 
channels respectively. The vertical grey lines are the corrected integer conductance 
values (see text), (b) Conductance histogram obtained from the same data set. The 
peak in the conductance histogram at Go extends to 53000 on the y-scede. (Inset) 
Schematic diagram of the configuration used in the analysis. The dark lines with 
arrows show the paths, which contribute to the conductance fluctuations in lowest 
order. (From Q) 



fluctuations, and from this description it is clear that the fluctuations are 
expected to vanish either when t n = 0, or when r n = 0. 
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Elaborating this model Ludoph et al. obtained the following analytical 
expression for acv , 



hkF vf \/l — cos 7 






a GV = ^ 2 - 71 l% =r^i) />-T2(l-T n ), ,:,;, 



where &f and «f are the Fermi wave vector and Fermi velocity, respectively, 
Tc — le/vF is the scattering time. The shape of the contact is taken into 
account in the form of the opening angle 7 (see the inset in Fig. [17|), and 



T^n is the applied voltage modulation amplitude. The full curves in Fig. [1Tb, 
are obtained from (g|), assuming a single partially-open channel at any point, 
i.e., assuming that channels open one-by-one as the conductance increases. In 
agreement with the results discussed above, the conductance for the smallest 
gold contacts is very well described by this simple approximation. The ampli- 
tude of the curves is adjusted to fit the data, from which a value for the mean 
free path is obtained, l c = 5 ± 1 nm. Similar experiments p2|J74| for copper 
and silver and for sodium also show the quantum suppression of conductance 
fluctuations observed here for gold, while for aluminium or niobium it is not 
observed. 

The back-scattering, which produces the conductance fluctuations, has a 
second observable effect, namely on the average conductance, as opposed to 
the fluctuating part. This effect is seen as a shift to lower conductance of the 
minima in agy, and can be described to lowest order by an effective resistance 
in series with the contact. This shift has been taken into account in the 
calculation of the curves in Fig. [l7|a. Also in the interpretation of conductance 
histograms, a phenomenological series resistance is often taken into account 
in order to describe the shift of the peaks to lower values p3,OJ36lp^] . The 
estimate of the mean free path obtained from this series resistance is indeed 
consistent with the value obtained by fitting expression (0) to the data. 

Thermopower. Further information on the quantum transport properties 
of atomic size contacts, again for gold, has been obtained from thermopower 
experiments using a modified MCB technique [J75|. By applying a constant 
temperature difference over the contacts, the thermally induced potential 
could be measured simultaneously with the conductance. Large thermopower 
values were obtained, which jump to new values simultaneously with the 
jumps in the conductance. The values are randomly distributed around zero 
with a roughly bell-shaped distribution, in contrast to what was expected 
from elementary free electron gas models |7q |. The thermopower was shown 
to result from the same mechanism as the conductance fluctuations, and 
the theory was adapted to describe the results. The experimental results 
were found to follow the law obtained from this defect-scattering model, and 
quantum suppression of the thermopower was indeed observed. By fitting 
the curves, a value for the mean free path was obtained which was in close 
agreement with the value from the conductance fluctuations. Since the two 
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experimental techniques are very different, and the typical energy scales for 
the effects are at least an order of magnitude apart (this scale is set by 
the modulation voltage amplitude, 20 mV, in one case and the temperature, 
~10 K, equivalent to ~ 1 mV, in the other), these results give strong support 
for the description and interpretation presented above. 

6.4 Implications for Conductance Curves and Histograms 

Let us first focus on results for the monovalent metal, gold, which has been 
most studied. The results obtained from superconducting subgap structure, 
shot noise measurements, conductance fluctuations and thermopower are in 
agreement and show that the conductance for the smallest contacts of mono- 
valent metals (gold) is carried by a single mode. For increasing contact size 
it is found that the transmission for the first mode goes to unity, before the 
second mode opens, then the second mode goes fully open, before the third 
opens, and so on. This property of atomic size contacts has been described 
as "saturation of the channel transmission" |72J . This property holds to very 
good approximation (a few percent) up to G — IGo, but deviations increase 
to 20% admixture of the next channels at G = 4Go. 

It is interesting to compare the positions of the maxima in the conduc- 
tance histogram and those for the minima in (Jgv in Fig. \\J\. It appears that 
these positions do not all coincide, which is most evident for the peak in 
the histogram at about G = 1.8Go- The histograms give preferential conduc- 
tance values, which may reflect a quantisation effect in the conductance as 
a function of contact diameter, but also a preference for forming contacts of 
certain effective diameters. Such preferential contact diameters may be ex- 
pected based on the fact that the contact is only a few atoms in cross section, 
which limits the freedom for choosing the diameter. It appears that at least 
the peak at 1.8 Go in the histogram for gold arises from this atomic geome- 
try effect. Although the shot noise and conductance fluctuation experiments 
both show that the conductance for gold contacts with G ~ 2Go is carried by 
two nearly perfectly transmitted modes, this conductance is not preferred, as 
evidenced by the conductance histogram. 



This observation amplifies the arguments given in the discussion of Fig. 11 
where it was argued that it is not appropriate to describe the conductance 
for gold as being quantised. There is no pronounced preference for conduc- 
tances near multiples of Go, with the exception of 1 Go, which is due to a 
special property of gold contact in that it forms chains of atoms at the last 
stages of contact elongation (see below) . Of course, the conductance is a true 
quantum property in the sense that it is carried by only a few well-defined 
modes, but the transmission for these modes can add up to any value of the 
total conductance. Moreover, we identify a property for these contacts, called 
the saturation of channel transmission, which describes the tendency for the 
conductance modes to open one-by-one. 
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For s, s-p and s-d metals the results for the subgap structure on the 
smallest contacts agree with the predictions of the tight binding model of 
Cuevas et al. under the assumption that the last plateau, before the jump 
to tunnelling, consists of a single atom contact. The agreement between the 
predicted number of channels and the number obtained from the fits of the 
current -voltage curves, thus confirm the assumption that the contacts consist 
of a single atom at the last stage of the contact elongation. The results on shot 
noise and conductance fluctuations for Al, where no quantum suppression was 
found, confirm the results by Scheer et al. that the conductance near G = 1 Go 
cannot be described by a single mode. We conclude that a saturation of 
channel transmission is not observed for Al, and is expected to be absent in 
all metals other than simple s-metals. Despite this apparent lack of a simple 
quantum mode structure, the conductance histogram for aluminium shows 
pronounced peaks near the first three or four multiples of Go |Q, albeit 
considerably shifted from perfect integers. A natural interpretation for the 
first peak in the histogram, which is consistent with the findings above, would 
be that it arises from a reproducible last contact configuration of a single 
atom, with a conductance close to 1 Go and three channels involved in the 
conductance. In this sense, the histogram peak would arise purely as a result 
of atomic structure. In analogy, the other peaks are also expected to relate 
to atomic structure, although it is less evident what this structure should be. 
In the simplest approach one could imagine that the peaks correspond to 1, 
2, 3 and 4 atoms in the contact cross section. 

Further support for the model of valence orbitals as a basis for the con- 
ductance channels comes from the anomalous slope of the plateaux of the 
conductance for aluminium, where conductance increases as the contact is 
stretched ( Fig. Q). Calculations of the conductance through a single atom, 
within the tight binding model discussed above, as a function of the bond dis- 
tance of the atom with its neighbours, reproduce the anomalous dependence 
of conductance on distance (77]]. For gold a nearly flat dependence is found, 
while for Pb the conductance decreases with increasing elongation, again in 
agreement with the observations. The anomalous slopes for Al are also found 
for larger contacts, and have been explained in terms of the stress-dependence 
of the electronic band structure [1f8| . It appears that this property of the band 
structure is conserved in contacts down to the atomic scale. 



7 Chains of Atoms 

All evidence shows that for a single atom contact for monovalent metals the 
current is carried by a single mode, with a transmission probability close 
to one. Guided by this knowledge, in experiments on gold Yanson et al. 
[ |79| discovered that during the contact breaking process the atoms in the 
contact form stable chains of single atoms, up to 7 atoms long. Independently, 
Ohnishi et al. |8(J] discovered the formation of chains of gold atoms at room 
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Fig. 18. Atomic configuration obtained during the last stages of breaking of a 
nanowire in a molecular dynamics simulation for gold, using a bath temperature of 
12 K. (From @) 



temperature in a combined STM and transmission electron microscope, where 
an atomic strand could be directly seen in the images. Currently, the only 
material for which this chain formation is observed is gold. It is interesting 
that for silver the effect is much less pronounced and copper does not seem 
to show it at all. 

Some understanding of the underlying mechanism can be obtained from 
molecular dynamics simulations. Already before the experimental observa- 
tions, several groups had observed the spontaneous formation of chains of 
atoms in computer simulations of contact breaking ]8l| , |82| , |33] ■ Figure O shows 
the results obtained by S0rensen et al. for gold. The authors caution that the 
interatomic potentials used in the simulation may not be reliable for this un- 
usual configuration. However, from these studies we may learn more about 
the atomic configuration sequences producing the chains, and they form a 
good starting point for more advanced model calculations. 

Such chains constitute the ultimate one-dimensional metallic nanowires. 
The current is carried by a single mode, with a transmission probability, 
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Conductance [in units 2e /h] 

Fig. 19. Histogram of the number of times each conductance is observed versus the 
conductance in units of the conductance quantum, Go, for sodium at T=80 K and 
bias voltage V=100 mV, constructed from over 10,000 individual scans. The smooth 
background (dashed curve) helps to bring out the smaller amplitude oscillations. 
(From ||) 



which is somewhat below 1 due to back-scattering, but can be tuned to unity 
by adjusting the stress on the junction. The chains sustain enormous currents, 
up to 80 |iA, due to the ballistic nature of the electron transport J7{|. More 
work is needed to elucidate the mechanism of chain formation, what limits the 
length of the chains, and why it works best for gold. Further unusual atomic 
configurations may be found, according to model calculations by Giilseren 
et al. Q, who show a series of "weird wire" structures depending on the 
number of atoms in the cross section of the wires and on the metallic element 
involved. 

8 Quantum Forces and Shell Structure in Alkali 
Nanowires 



Apart from their role in determining the electronic transport properties of 
the contacts, the quantisation of the wave functions may also affect the en- 
ergy of formation of the contact. In direct analogy with the effect of elec- 
tronic shell closing on the formation energy of clusters |31|] , one expects that 
specific contact diameters will be stabilised by the formation of quantum 
modes. Several groups have recently presented model calculations on this 
problem, mostly considering free electrons confined by a hard wall potential 
[ p6p7| , |88| , ^9U9C|j9l[ | . Local density calculations taking the atomic structure 
into account have also appeared |92,p3,p3. There appears to be a consensus 
on the magnitude of the force fluctuations which may result from this mech- 
anism, which is of the order of 1 nN, and is comparable to the force jumps 
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for the smallest contacts observed in Fig. 0. Some authors argue that the 
jumps observed in the conductance as a function of contact stretching should 
therefore be considered as being the result of the underlying electronic quan- 
tum modes. However, in the experiment many conductance steps are found 
which are much smaller, or much larger, than a conductance quantum, and 
all steps show similar mechanical and dynamic behaviour. In addition, the 
conductance steps in monovalent metals look similar to those in sp or sd 
metals, where such simple description definitely breaks down. Experimental 
methods will have to be developed to test a possible quantum-mode-based 
atomic force mechanism. One possibility is to measure the charge, or the 
work function of the nanowire, which is expected to fluctuate in unison with 
the force within this simple free electron gas picture |87||KJ • 

While the contribution of quantum modes to the force in the smallest 
contacts is still under debate, a newly observed phenomenon for larger con- 
tacts provides strong evidence for quantum force fluctuations. In a recent 
study of conductance histograms for sodium, potassium and lithium up to 
conductances much larger than shown in Fig. 0, Yanson et al. p3] observed a 
large number of additional peaks (Fig. p9|) . The peaks are not as sharp as the 
ones associated with conductance quantisation at low conductance (Fig. ||) 
and cannot be identified with multiples of the conductance quantum. The 
peaks become more pronounced as the temperature is raised to about 80 K, 
and the position of the peaks is seen to be periodic in the square root of the 



conductance, as illustrated in Fig. 20 



The interpretation of the phenomenon is based on fluctuations in the 
density of states as a function of energy (or equivalently as a function of 
diameter) for a free electron gas inside a cylindrically symmetric wire. At the 
points where a new mode (the bottom of a 1-dimensional subband) crosses 
the Fermi energy, the density of states shows a l/\/E — E n singularity. These 
singularities are smeared out by the finite length of the wire and would not 
have a very pronounced effect were they homogeneously distributed. However, 
the symmetry of the wire gives rise to a bunching of the singularities, which 
can be associated with the electronic shells in metal clusters. The resulting 
density of states fluctuations have been analysed by Stafford and coworkers 
[ p6p0|| , by Yannouleas et al. |58| and by Hoppler and Zwerger |pi[ . The 
fluctuations in the density of states result in local minima in the total energy 
for the nanowire, and the stable wire diameters predicted from this model 
are in fairly good agreement with the observed periodic peak structure in the 
histograms f85fl . 

It was found that a direct comparison with cluster magic numbers is pos- 
sible. One can calculate the effective cluster radius from the magic numbers 
for sodium clusters ||l[] N m = 8, 20, ..., 1500 using R m — rgiV m , where r$ 
is the Wigner-Seitz radius of the atom. Then the conductance of a nanowire 
with these magic radii can be found using the semi-classical expression p3| , 
G m — Go(TrR m /\-p) 2 (l — Af/7t_R to ). The square root of the conductance 
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Fig. 20. Square root of the positions of the maxima in Fig. Il9|, yjGi/Go, versus 
the shell number {open squares). The experimental peak positions are compared to 
the expected positions derived from the magic numbers of metallic clusters [closed 
circles). (From |8q]) 



values obtained in this way are plotted against the shell number m and com- 
pared to the histogram peak positions in Fig. KG. A striking agreement is 
observed, which is believed to result from the fact that the dominant fluctu- 
ation terms in the density of states can be described in terms of the lowest 
order semi-classical trajectories inside the system p5| , p6||9l|] . These trajecto- 
ries (diametrical, triangular and square orbits inscribed inside the sphere and 
cylinder, respectively) are the same for clusters and nanowires. However, the 
agreement in Fig. E(J is better than expected since the relative contribution 
of each of the trajectories for the two systems should be different. Further 
work is needed to clarify the details of the shell structure in nanowires. 



9 Discussion and Outlook 

With the developments of the last few years a microscopic understanding of 
atomic-scale electrical transport properties is beginning to grow. A coherent 
picture of conductance modes in a single atom derived from atomic valence 
orbitals is obtained, which finds strong support in the experimental observa- 
tions for various metallic elements. For contacts of several atoms, monovalent 
metals demonstrate saturation of channel transmission, which is a very in- 
teresting experimental observation, but is not yet fully understood on the 
microscopic level. Also, the nature of peaks in the histograms is still not fully 
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resolved: quantisation plays a role but not exclusively. This is evident from 
the histograms for aluminium, in combination with the channel numbers ob- 
tained from subgap structure, shot noise and conductance fluctuations. The 
peak at 5 Go in sodium is stronger than expected from a simple free electron 
nanowire model, and the peak below 2 Go in gold does not correspond to 
the position of the conductance quantum, even after correction for a series 
resistance, as was seen in conductance fluctuation experiments. 

It seems clear that a full understanding of the structure in conductance 
histograms requires a description of the mechanical evolution of the atomic 
structure of the contacts. It is possible that conditions such as temperature or 
bias voltage must be taken into account |97J . Much has been learned already 
from molecular dynamics calculations, but new experimental tools are prob- 
ably essential. Promising developments come from two groups, which have 
recently shown transmission electron microscopy images with atomic resolu- 
tion of contacts that can be controlled in situ |80J . |98| . p9[ . By simultaneous 
measurement of the conductance while observing a single atom contact in 
the electron micrograph, Ohnishi et al. clearly confirmed that a single atom 
gold contact has a conductance close to 1 Go. 

The question of the effect of the conductance modes on the force in metal- 
lic contacts will, without doubt, receive a lot of attention in the near future. 
A related question at present being considered is whether the increased den- 
sity of states at specific wire diameters may be lifted by deformation of the 
shape of the wire cross section (i n ana logy to the Jahn-Teller effect) |87J] or 



by a spontaneous magnetisation 1 100 . Many other open questions promise 



new and interesting developments in this area of physics. The dynamics of 
atomic structures, seen as two-level fluctuations, have not received sufficient 
attention. A related problem concerns the interaction with the electric cur- 



rent and the heating of the contact [ 101 1 . The observed chain formation for 



gold may lead to a deeper understanding of bonding forces on the atomic 
scale. We may be able to produce much longer atomic wires, which will open 
a new field for study of purely one-dimensional solid state physics. 
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